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Seismic  Source  Depths  from  Depth  Phase^. 
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The  primary  concept  involved  in  the  improved  techniques  is  the  util- 
ization of  source  depth  information,  contained  throughout  the  seismi 
coda,  which  has  previously  either  been  ignored  or  improperly  used. 
This  information  was  utilized  through  the  computation  of  travel 
times  for  t lie  later  arriving  seismic  phases  and  the  development  of 
techniques  which  improve  the  detection  of  these  later  phases. 

The  results  of  this  research  demonstrated  that  the  automated  source 
depth  determination  procedure  constructively  utilizes  source  depth 
information  contained  throughout  the  seismograms,  and  thereby  sig- 
nificantly enhances  source  depth  determinations.  This  wa s shown  by 
obtaining  the  correct  depth  for  an  event  through  the  separate  analy- 
sis of  the  first  and  second  minute  of  data,  verifying  the  usefulness 
of  source  depth  information  contained  in  the  later  seismic  phases. 

A further  verification  of  the  effectiveness  of  this  procedure  was 
achieved  by  obtaining  correct  depth  estimates  from  data  arriving 
through  individual  seismic  phases,  and  from  data  recorded  at  indi- 
vidual stations. <- 


During  this  contract,  research  was  also  coordinated  which  establish 
ed  possible  modifications  to  these  techniques  which  should  enhance 
their  effectiveness  for  events  having  shallow  source  depths. 


SeCOHIty  CL  ASStOCATIOM  OF  This  t'Mit  OFh.n  »«l*  Knr«r»rf) 


I 


SUBJECT: 


Improved  Procedures  for  Determining  Seismic 
Source  Depths  from  Depth  Phase  Information 


t 

AFTAC  Project  No 

ARPA  Order  No 

ARPA  Program  Code  No 

Name  of  Contractor 

Contract  Number 

Effective  Date  of  Contract 
Report  Period 

Amount  of  Contract 

Contract  Expiration  Date.. 
Project  Scientist 

t 


VELA  T/6710 

2551 

6F1 0 

ENSCO,  INC. 

F08606-76-C-0003 

September  1,  1975 

September  1,  1975 
to  June  30,  1976 

$59,998 

June  30,  1976 

Edward  A.  Page 
(703)321-9000 


t 


I 


TABLE  OE  CONTENTS 


Page 

1.0  INTRODUCTION  1-1 

2.0  AUTOMATED  SEISMIC  SOURCE  DEPTH  ANALYSIS  2-1 

2.1  Seismic  Source  Depth  Analysis  Procedures  2-1 

2.2  Access  and  Storage  of  Differential  Travel  Time 

Inform  ition  2-5 

3.0  APPLICATION  OF  THE  ANALYSIS  PROCEDURE  TO  THE  ILLINOIS 

EVENT  3-1 

4.0  SHALLOW  SHOURCE  DEPTH  DETERMINATION  USING  CEPSTRUM 

ANALYSIS  4-1 

4.1  Analysis  of  the  Boxcar  Event  4-1 

4.2  Determination  of  the  Effectiveness  of  Cepstrum 

Analysis  for  Shallow  Source  Depths  Using 
Synthesized  Data  4-10 

5.0  CONCLUSION  5-1 

APPENDIX  A:  COMPUTATION  OF  THE  CEPSTRUM  A-l 

APPENDIX  B:  COMPUTATION  OF  THE  CEPSTRUM  MATCHED  FILTER  B-l 


» 


1.0  INTRODUCTION 


The 

ob  j 

ect i ve 

o f 

this 

research  wa 

s to 

formul 

at 

e se i 

ism 

i c 

SOU! 

:cc  de 

pth 

deter 

min  a 

tion 

teqhniques , 

deve 

loped 

du 

ring 

th 

i s 

and 

two  p 

rev 

i ous  c 

ontracts  , 

into  an  au 

tomated  ana 

lv 

sis  procedu 

re 

The 

purpo 

se 

is  to 

subs 

tanti 

ally  increa 

se  th 

e perc 

en 

tage 

of 

eve 

nt 

for 

which 

ac 

curate 

sou 

rce  depths  can  be  dot 

ermine 

d 

and  1 

*o 

faci 

li 

tate  the 

ana 

lysis . 

In 

this 

procedure , 

the 

source 

d 

epth 

is 

det 

er 

mined  by  the  degree  to  which  cepstrum  patterns,  computed  from 
different  portions  of  single  or  multistation  seismic  data,  agree 
with  those  patterns  expected  for  a given  source  depth.  By  auto- 
matically accounting  for  variations  in  the  differential  travel 
times  caused  by  different  station  locations  and  the  presence 
of  later  propagation  modes,  the  depth  phase  information  contained 
throughout  the  seismograms  contribute  to  this  depth  estimate. 


During  this  contract,  the  logic  and  algorithms  necessary  to 
incorporate  and  automate  the  various  source  depth  analysis  techn; - 
ques  were  developed  such  that  the  seismic  data  and  station  to 
source  distances  are  specified  to  obtain  this  source  depth  estimate. 
This  work  included  the  development  of  algorithms  needed  to  store 
and  access  differentia!  travel  times  as  a function  of  source  depth 
and  station  to  source  distances,  for  several  different  propagation 
modes  . 


The  effectiveness  of  this  automated  source  depth  determination 
procedure  was  demonstrated  through  its  application  to  the  Illinois 
Event.  The  results  of  this  analysis  showed  that  the  depth  phase 
information,  contained  throughout  the  seismograms,  was  in  fact 
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being  constructively  utilized  and  thereby  enhanced  the  depth 
estimate.  This  was  shown  by  obtaining  the  correct  depth  for 
this  event  through  the  separate  analysis  of  the  first  and  second 
minute  of  data  verifying  the  usefulness  of  depth  phase  informa- 
tion contained  in  the  later  seismic  propagation  modes.  A further 
verification  was  achieved  by  obtaining  this  depth  estimate 
through  the  analysis  of  data  arriving  through  individual  propa- 
gation modes.  In  addition,  it  was  shown  that  accurate  depth 
estimates  could  be  obtained  from  the  analysis  of  data  recorded 
at  single  stations.  The  primary  result  of  these  analyses  is  that 
unlike  other  analysis  procedures,  this  automated  seismic  source 
depth  determination  procedure  is  utilizing  depth  phase  information 
contained  throughout  the  seismograms  and  should  allow  one  to 
obtain  accurate  depth  estimates  for  a higher  percentage  of  events 
then  previously  possible. 

During  this  contract,  research  was  also  conducted  to  investi- 
gate the  effectiveness  of  these  techniques  in  determining  depths 
of  shallow  event;  and  to  determine  what  modifications  would  enhance 
the  analysis  for  this  type  of  data.  This  research  involved  the 
analysis  of  both  the  Boxcar  Event  and  synthesized  shallow  depth 
data . 
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2.0  AUTOMATED  SEISMIC  SOURCE  DEPTH  ANALYSIS 


In  order  to  both  enhance  and  automate  seismic  source  depth 
determinations,  various  techniques  developed  during  this  and  two 
previous  contracts,  were  formulated  into  a computerized  procedure 
allowing  the  analyst  to  input  seismic  data  and  station  to  source 
epicenter  distances  to  obtain  a source  depth  estimate.  In  this 
procedure,  the  source  depth  is  determined  by  the  degree  to  which 
cepstrum  patterns,  computed  from  different  portions  of  multi  or 
single  station  seismic  data,  agree  with  those  expected  for  a given 
source  depth.  The  variations  in  differential  travel  time  caused 
by  the  different  station  to  source  distances  and  the  presence  of 
later  propagation  modes  are  both  accounted  for  and  allow  depth 
phase  information  contained  throughout  the  seismogram  to  contribute 
to  this  depth  estimate. 
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2.1  Seismic  Source  Depth  Analysis  Procedure 

A prototype  of  the  seismic  source  depth  analysis  procedure  is 
flow  charted  in  Figure  2.1.  The  input  data  consists  of  the  seismic 
data  recorded  over  a suite  of  stations,  the  data  sample  length 
(governing  the  maximum  differential  delay  time  observable),  the 
length  of  coda  to  be  analyzed,  the  offset  times  between  the  start 
of  the  data  and  the  P-wave  onset,  and  the  station  to  source  epicenter 
distances . 


The  automated  analysis  then  proceeds  to  select  a data  sample 
from  a given  station  recording  and  coda  position  and  a cepstrum 
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and  ccpstrum  matched  filter  (CMP)  output  is  then  computed  from 
this  data  sample  (Appendix  A and  Appendix  B give  details  of  the 
steps  in  computing  the  ccpstrum  and  CMP).  A flow  chart  of  the 
CMP  is  shown  in  Figure  2.2.  The  first  of  a range  of  trial  source 
depths  is  then  selected.  For  this  trial  depth,  and  station  to 
source  epicenter  distance,  the  differential  travel  time  for  the 
following  propagation  modes  are  accessed  from  computer  storage: 
pP-P,  PP-P,  pPP-PP,  PPP-P,  pPPP-PPP,  PcP-P,  and  pPcP-PcP.  Prom 
the  start  and  end  time  of  the  data  sample,  it  is  then  decided 
which  propagation  modes  will  contribute  to  the  CMP  output  for  this 
particular  data  sample,  and  at  what  surface  reflected  delay  times. 
The  maximum  CMP  output  within  a given  tine  window  of  this  expected 
delay  time  is  accumulated  separately  for  each  contributing  pro- 
pagation mode.  For  this  same  data  sample,  the  procedure  is  repeated 
for  each  of  the  trial  depths  to  be  assumed.  This  completes  the 
analysis  of  this  data  sample  and  the  next  sample  is  likewise 
processed.  For  each  trial  depth  the  CMF  output  is  accumulated 
over  the  different  data  samples  for  each  propagation  mode.  The 
accumulated  CMP  outputs  for  each  propagation  mode  are  then  scaled 
and  summed  such  that  each  mode  has  equal  contribution  and  the 
final  result  is  a plot  of  the  cumulative  CMP  output  versus  depth. 

The  cumulative  CMF  output  versus  depth  for  each  propagation  mode 
is  also  available  as  an  optional  output. 
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2 . 2 Access  and  Storage  of  Differential  Travel  Time  I n forma  t i on 

In  order  to  perform  the  automated  source  depth  analysis  one 
needs  to  determine  which  seismic  propagation  modes  have  arrived 
during  a given  data  sample  and  what  are  the  differential  travel 
times  between  these  arrivals  and  their  associated  surface  reflec- 
tions. For  the  event  distances  and  depths  of  interest  in  this 
work  the  following  seismic  propagation  modes  are  involved:  pP-P, 

PP-P,  pPP-PP,  pPPP-PPP,  PPP-P,  pPcP-PcP  and  PcP-P.  The  automated 
analysis  needs  to  access  differential  travel  times  for  each  of 
these  propagation  modes  over  the  range  of  depths  from  0 to  100  km 
and  over  epicenter  distances  of  10°  to  90°  to  an  accuracy  of  a 
few  tenths  of  a second  in  most  cases. 

The  travel  time  differences  were  obtained  by  the  application 
of  ray  tracing  using  the  spherically  symnetic  isotropic  earth 
velocity  model  used  for  the  BSSA  seismological  tables.  A three 
dimensional  polynominal  surface  representation  of  the  differential 
travel  time  was  used  to  facilitate  computer  access,  reduce  computer 
storage  load  and  to  perform  the  necessary  interpolation  of  the 
computed  values.  We  now  describe  the  procedure  used  for  obtaining 
the  representation  of  the  differential  travel  times  as  a function 
of  source  depth  and  source  to  receiver  distance  for  the  following 
seismic  propagation  modes:  pP-P,  PP-P,  pPP-PP,  PPP-P,  pPPP-PPP, 

PcP-P  and  pPcP-PcP. 
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Case  I:  Surface  Fits  to  pP-P,  pPPP-PPP,  pPP-PP,  pPcP-PcP 


An  examination  of  Figure  2.1.1  suggests  we  represent  these  data 
as  a double  power  series  in  which  the  travel  time  difference  t 
is  the  dependent  variable  and  the  depth  d and  epicenter  distances 
A are  the  independent  variables.  The  surface  of  travel  time 
differences  is  then  given  by 


l, A)  = 7 T.  A^d1 

iJ 


The  coefficients  t..  and  the  number  of  them  needed  were  determined 

iJ 

as  follows.  Each  curve  at  constant  Ac  was  found  to  be  adequately 
represented  by  a cubic  in  d: 

T(d,Ac)  = T1(Ac)d  + T2(A(.)d2  + T3(Ac)d3 


Here  t vanishes  at  zero  depth.  The  least  squares  values  of  t.(a) 
were  then  in  turn  adjusted  to  a power  series  in  A: 

N 

T i ( A ) = £ t A^  , i = 1 » 2 , 3 

j=o 

Reliable  input  data  for  A less  than  10°  were  not  available,  so  no 
attempt  was  made  to  force  x ^ ( A ) to  vanish  at  zero  epicenter  distance 
Thus,  our  representation  should  be  used  only  for  A >y  10°.  All  the 
data  sets  of  Case  I can  be  represented  by  N=9.  Examples  of  the 
effectiveness  of  this  representation  are  shown  in  Figures  2.1. -2a 
through  2.1. 2d  for  pP-P,  pPP-PP,  pPPP-PPP  and  pPcP-PcP. 
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Plot  of  travel  time  differ 
ence  (seconds)  for  pP-P 
vs.  depth  (kni)  for  given 
epicenter  distances  (A0) 


Depth  (km) 

Figure  2.1.1 
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Figure  2.1.2a.  Differential  Travel  Times  Obtained  from  Polynominal  Curve  Fit  for  nP-P 


I 


O 

•AJ 

a jo 

vf 

4 

r-* 

in 

O 

a 

LJ 

x. 

X 

rvj 

Oj 

OJ 

o 

P- 

p-  in 

rv 

X 

CP 

o 

> 

O 

X 

—4 

X 

PJ 

*a4 

OJ 

n 

-n 

<1 

o 

4 

X 

•l 

Na 

7 

o 

o 

O 

z> 

■Zi 

o 

G 

^ ■ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

>oJ 

• 

• 

« 

• 

• 

• 

• 

• 

• 

1 

o 

^-4 

^y 

X 1 

X 

o 

A.' 

p- 

ro 

o 

o 

1 

1 

1 

i 

1 

jJ 

CP 

r-4 

—4 

»— 4 

AJ 

2 

7 

3 

UJ 

1 

i 

a 

2 

o 

vX 

X 

y 

n 

X 

P- 

•4- 

UJ 

o 

X 

y 

ro 

7 

— H 

ro 

a 

* - 4 

n 

a 

in 

X 

• 4 

D 

in 

X 

UJ 

o 

•wO 

a x a 

a 

— « 

ro 

o 

LU 

o 

r- 

r- 

OJ 

o 

—* 

3, 

'G 

o 

o 

o 

o 

o 

o 

*-4 

—4 

r^4 

■> 

o 

XI 

-y 

x - 

ro 

X 

X 

X 

p- 

Lu 

o 

o 

o 

13 

a 

o 

o 

—4 

—a 

1 

1 

i 

i 

1 

1 

1 

1 

1 

1 

• 

• 

• 

• 

• 

• 

• 

• 

• 

•—4 

• 

• 

• 

• 

• 

• 

• 

• 

• 

o 

UJ 

a1 

aJ 

a 

su! 

uJ 

UJ 

aJ 

UJ 

Jj 

o 

o 

^4 

4- 

X 

X 

o 

AJ 

o 

AJ 

a 

o 

i 

I 

1 

1 

1 

LU 

3 

y 

X 

— < 

N- 

p- 

CP 

X 

X 

a4 

-H 

rH 

a 

X) 

O 

ro 

a 

o 

X 

p^ 

X 

X 

•PJ 

X 

ro 

a 

UJ 

OC 

a 

n 

p- 

x 

X 

X 

s0 

o 

2 

LU 

a 

x 

*o 

O 

X 

»a 

n 

■J 

X 

X 

o 

LU 

o 

O' 

X 

4*  - 

r— 4 

p- 

X 

o 

P- 

►“4 

o 

X 

rv 

f— 4 

a 

o 

y 

o 

PJ 

< 

U 

3 

■jr 

X 

sC’ 

rp 

X 

x 

o 

O' 

ro 

o 

•m 

—4 

X 1 

AJ 

X 

r*"s 

X 

V— 

o 

p- 

pj 

r— 4 

n 

Pa 

o 

ro 

ro 

2 

CO 

X 

X 

X 

o 

p~ 

a* 

n 

"0 

in 

>■■■4 

~r~ 

o 

<y 

X 

PJ 

lO 

X 

in 

r-H 

1 

o 

o 

o 

r— 4 

zz 

— 

o 

—4 

o 

3 

UJ 

pa 

a 

cp 

X 

ro 

O 

O' 

o 

m 

U- 

• 

• 

• 

• 

• 

• 

• 

• 

• 

Li 

• 

• 

• 

' • 

• 

• 

• 

• 

• 

U' 

x 

3 

pa 

l\l 

'J 

■O 

a 

pp 

'XJ 

>3- 

a 

o 

»-4 

•y 

n 

CG 

o 

AJ 

vn 

a 

UJ 

o 

1 

i 

1 

i 

* 

1 

3 

UJ 

<3 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

uJ 

p> 

— * 

r— 

AJ 

> 

1 

u_ 

1— 

» 

1 

1 

i 

1 

3 

< 

Lu 

u. 

1 

a 

o 

; 

Q 

o 

m 

—4 

p-  « 

a 

»— 

O 

m. 

X 

h- 

o 

—4 

X 

X 

ro 

o 

r 

2 

c 

! 

LU 

o 

3 

X 

X X 

X 

7 

X 

1 

o 

rc 

o 

ro 

o 

o 

o 

o 

<1 

o 

S) 

ro 

y i 

r-*4 

n 

4 

O' 

Q 

o 

o 

o 

O 

o 

n 

«.*o 

— « 

—4 

1 

UJ 

O 

• 

• 

• 

• 

• 

• 

• 

• 

• 

UJ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

X 

21 

3 

o 

^-4 

>y 

n x 

o 

X- 

►— 

o 

i 

1 

1 

1 

1 

1 

»~4 

fVJ 

in 

X 

O' 

—4 

t 

£ 

3 

3 

—4 

r* 4 

AJ 

<1 

X 

o 

1 — 

o 

o 

o 

o 

g 

O 

—a 

—4 

<1 

1 

1 

1 

» 

I 

1 

i 

1 

1 

1 

1 

2 

3 

M 

_i 

U1 

uj 

! • ! 

UJ 

UJ 

UJ 

UJ 

UJ 

LJ 

uJ 

r> 

o 

CG 

'O  T'  i 

ro 

ro 

ro 

O 

O 

o 

P~ 

rvj 

•5 

X 

Pa 

p- 

7 

Va1 

0. 

UJ 

3 

—4 

—4 

c 

p- 

4- 

o 

X 

X 

ro 

*“4 

o 

O'  x • 

—4 

P- 

7 

OJ 

AJ 

3 

o 

p- 

r— 4 

7 

y 

o 

rv 

o 

•n 

Lul 

> 

a 

(\J 

'n 

O ; 

p- 

o 

> 

o 

-4 

ro 

Pa 

cn 

o 

n 

rv 

ro  < 

o 

r-» 

r- 4 

X 

7 

< 

o 

o 

o 

o 

G 

OJ 

G 

o 

OJ 

<1 

t/ 

IT 

r 

a 

a 

1/ 

fV 

c 

LT 

X- 

• 

• 

• 

• 

• 

• 

• 

• 

• 

C' 

• 

• 

• 

• 

• 

• 

• 

• 

• 

c 

X 

a 

PJ 

X* 

X 

3 

ro 

o 

T 

ro 

< 

o 

»-4 

•y 

X 1 

X 

o 

As. 

in 

O 

o 

1 

1 

1 

i 

1 

2 

h- 

3 

X 

V 

o 

CO 

o 

* 

O' 

n 

lO 

X 

u' 

m 

* 4 

rH 

r— 4 

AJ 

3 

in. 

4 

oO 

'"O 

X 

pa 

X 

Ij' 

-*-4 

j' 

_rv 

J 

2 

o 

j 

X' 

o 

<y 

■c 

• — 

o 

CP 

o 

1 

<1 

2 

Pa 

ro 

r 

r— 4 

(V 

rw 

y 

<« 

o 

C 

P-  7 

c 

A* 

»— 4 

X 

0 

o 

fa. 

Os 

1 

Pa 

— • 

• — 

•— 1 

G 

31 

»~4 

<1 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

4* 

o 

o 

a x 

IP 

-o 

—4 

X 

X 

K- 

o 

— 

X 

—4 

o 

O 

J 

C 

1 — 

i 

1 

i 

1 

i 

1 

o 

y 

O' 

p-  . 

r*4 

o 

r- 

X 

X 

3 

JO 

o 

o 

o 

o 

—4 

V 

<1 

* 

• 

• 

• 

• 

• 

• 

• 

• 

• 

a 

• 

• 

• 

• 

• 

• 

• 

• 

• 

3 

$ 

o 

*— 4 

ro 

X ! 

•a. 

CP 

o 

ro 

^1 

1 

1 

1 

1 

1 

T 

3 

<« 

4- 

»-4 

r»4 

r— 

U-4 

y 

C_> 

( 

1 

' — 

j 

a. 

-J 

UJ 

UJ 

< 

i 

UJ 

o 

-O 

CP  AJ  < 

o 

CP 

7 

A 

a 

O 

X 

n 

fa. 

X 

-o 

a 

T 

A 

o 

U 

^-4 

o 

» ’4 

OJ 

y 

■n 

O' 

r—4 

-y 

_J 

o 

r—*  (\1 

X 

y 

O 

A 

21 

o 

X 

n 

CP 

X 

i 

•JO 

X 

A 

o 

o 

o 

o 

o 

o 

O 

r-< 

X 

o 

fu 

o rv 

X 

r**4 

A 

lu' 

o 

C 

o 

o 

o 

PJ 

OJ 

A 

x 

a 

4 

4 

1 

» 

1 

1 

• 

i 

1 

9 

C 

• 

• 

• 

• 

• 

• 

• 

• 

A 

a 

• 

• 

• 

• 

• 

• 

• 

• 

A 

> 

o 

u. 

•u‘ 

Lu 

u. 

a- 

Lu 

UJ 

u: 

’'Uj 

li' 

K 

o 

1— 

ro 

i r 

■£ 

r- 

X 

o 

A 

u 

o 

i 

;« 

1 

c 

i 

A 

Ll. 

3 

X 

a 

f\ 

r— . 

P* 

X 

n 

rvj 

. — i 

7 

ro 

*» 

f— . 

A 

Lu 

m 

i 

A 

•— 

X 

*0 

T ' 

Pa 

a 

m 

r-» 

-o 

rvj 

rO 

•aJ 

A 

L- 

A 

o 

wf 

c 

0. 

P. 

•» 

*— 

X 

IT 

X 

4— 

•X 

U 

{ 

*—• 

i 

UJ 

I— 

X 

Pa 

CV 

O' 

a- 

LT 

X 

7s 

in 

2 ' 

o 

*-* 

a i 

X 

X 

s> 

A 

D 

o 

ro 

X 

cT* 

— 

^“4 

o 

A 

A 

X 

2 

3 

rp 

!> 

r- 

Pa 

X 

JO 

IT- 

Pa 

x 

UJ 

o 

O 

3 

X 

7 

•—* 

*7 

A 

o 

y 

in 

o 

y 

n. 

A 

A 

• 

UJ 

-0 

r' 

ro 

-r" 

• ■ 

O' 

ru 

m 

o 

o 

x> 

a ' 

— 

AJ 

x < 

I— 

m 

c* 

A 

UJ 

o 

O 

— 

r— 

Os« 

A 

A 

3 

X 

G 

rz 

n 

*.n 

a 

..p 

# ’» 

—4 

UJ 

• 

• 

• 

• 1 

i • 

• 

• 

-•> 

A 

2 

• 

• 

• 

' • 

• 

• 

• 

A 

A 

u. 

u 

3 

a 

r-a 

X) 

rvj 

rvi 

^4 

«\J 

u. 

“» 

*—4 

y 

X 

a 

X 

A 

H-4 

j> 

1 

1 

1 

A 

A 

u. 

< 

' • 

• 

• 

• 

• 

• 

• • 

• 

• 

0 

Lu 

AJ 

i 

A 

*- 

a 

1 

A 

A 

*-4 

u. 

u- 

l 

( 

I 

i 

1 

*— 4 . 

j 

A 

A 

A 

3 

a 

1 

3 i 

3 

Ul 

l 

1 

o 

O 

O X 

y 

•O' 

in 

‘4« 

A 

3 

O 

Pa 

r-4 

r— i 

rv 

PJ 

A 

A 

UJ 

o 

| 

i 

| 

UJ| 

o 

p^ 

-y 

y 

o 

X 

A 

> 

r> 

1 

y 

o 

O 

o 

o 

G 

A 

A 

X 

o 

i 

x • 

o 

3> 

^ X>  CO 

CZ* 

X 

A 

< 

o 

-V 

o 

v3 

X 

A 

A 

M 

i 

, 

*— • 

• 

• 

• 

• 

• 

• 

• 

A 

a 

• 

• 

• 

• 

• 

• 

• 

A 

A 

X 

! 

1 

i 

i- 

r4 

-o 

y x 

r- 

Ct 

A 

X 

A 

A 

lU 

1 

<—4 

A 

1 

•““4 

A 

A 

_j 

M 

3 

A 

77. 

A 

A 

UJ 

X 

o 

0~* 

•OJ  1 

3 

y 

n 

O 

X 

O' 

UJ 

lu) 

1 

> 

LU 

> 

o 

X 

>y 

A,  ' 

a 

X 

A 

U-. 

o 

' 

ro 

Pi 

o 

A 

A 

A 

< 

X 

<J 

o 

O' 

a x t 

n 

* 

A 

A 

*?: 

o 

X 

a 

a 

»n 

X 

A 

A 

A 

x 

11 

II 

II 

II 

II 

II 

II 

m 

II 

II 

a 

—4 

a 

n * 

•a 

3 

* 

A 

i— 

o 

o 

O 

o 

o 

o 

A 

A 

A 

►— 

rr 

V— 

• 

• 

• 

• 

• 

• 

$ 

A 

A 

(ul 

• 

• 

• 

• 

• 

• 

A 

A 

A 

UJ 

3 

3 

“)  ■ 

3 

3 

3 

3 

3 

^-4 

<0 

v4  . 

n 

P- 

* 

A 

A 

X 

o 

1 

A 

A 

A 

I 

3 

^-4 

A 

A 

— H 

A 

A 

A 

o 

Lu 

« 

4 

A 

in 

t 

A 

A 

A 

a 

a 

V—  . 

3 

i 

a 

< 

4 

-> 

J' 

r , 

r 

c* 

c_* 

J. 

<r 

t 

-* 

c. 

c. 

s~* 

- 

C 

ZT 

i 

UJ 

3 

o 

o 

c i 

rj> 

<7 

3 

w 

•u* 

,3 

"O 

v_3 

O 

a 

_l 

1 

<r 

T O 

O' 

O 1 

o 

O 

O 

O 

p 

< 

X 

o 

o 

•D 

rp 

*> 

O 

O 

a 

cr 

i 

* 

U 

• 

• 

• 

• 

• 

• 

• 

• 

»-4 

h- 

1“ - 

• 

• 

• 

• 

• 

• 

• 

• 

3 

3 

3 

a 

Ul 

*o 

y> 

O' 

o 

o 

o 

in 

3 

D. 

X 

0*1 

n 

"Zj 

o 

♦ 

O 

<1 

UJ 

3 

rv  r> 

"O 

■Jl 

n- 

o 

UJ 

UJ 

LU 

•—* 

a. 

ro 

•o 

X 

P- 

G 

a 

c 

1 

1 

1 

O 

c 

C 

i 

a 

c 

Q 

3^ 

J 

i 

1 

1 

i 

i 

l 

i 

» 

2-9 


Figure  2.1.2b.  Differential  Travel  Times  Obtained  from  Polvnomir.al  Curve  Fit  for  pPP- 
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Figure  2.1. 2d.  Differential  Travel  Tines  Obtained  fron  Pol vnomi nal  Curve  Fit  for  nPcP-PcP 


I 


Each  differential  travel  time  surface  considered  in  this 
work,  except  those  representing  core  reflections,  exhibit  an 
unallowed  region,  i.e.,  there  exists  a range  of  epicenter  distances 
for  which  below  a certain  depth  the  surface  reflected  mode  is  not 
received.  The  boundary  of  such  a region  projected  into  the  (d,A) 
plane  for  pP-P  is  shown  in  Figure  2.1.3. 


d (km) 


Figure  2,1.3 


These  boundaries  can  be  adequately  represented  by  a parabola  so 
that  our  polynomial  model  is  subject  to  the  constraint  that  for 
a given  A we  must  have  d < a+bA  , where  a and  b have  been  deter- 
mined from  least  square  adjustments  to  plots  like  that  of 
Figure  2.1.3.  In  Figure  2.1.2  the  unallowed  region  is  shown  as 
asterisk  table  entries. 

Case  II:  Line  Fits  to  PP-P,  PPP-P,  PcP-P 

These  travel  time  difference  plots  can  be  represented  as  single 
curves  depending  on  A only,  since  between  5 and  100  km,  the  depth 
dependence  is  very  weak.  Thus  we  have 


2-12 


The  PP-P  and  PPP-P  differences  appear  to  vanish  as  t\->o , so  for 
these  surves  we  set  t =0  and  also  find  N=8  is  sufficient.  Poi 
PcP-P,  the  core  reflection  insures  that  tqo  + 0.  Here  N=6  is 
sufficient . 

This  polynomial  surface  representation  of  the  differential 
travel  times  can  be  stored  using  a total  of  143  coefficients  and 
requires  less  than  .02  seconds  to  compute  these  times  for  eight 
propagation  modes  on  the  CDC  6600  computer. 


3.0  APPLICATION!  OP  THE  ANALYSIS  PROCEDURE  TO  THE  ILLINOIS  EVENT 


In  order  to  establish  that  this  automated  analysis  procedure 
was  in  fact  constructively  utilizing  depth  phase  information 
throughout  the  seismic  coda  anu  thereby  enhancing  one's  ability 
to  extract  accurate  source  depth  estimates,  the  Illinois  Event 
was  analyzed  in  several  ways.  Although  depth  phase  delay  times 
can  be  obtained  using  conventional  procedures  for  this  event 
by  demonstrating  that  these  delay  times  can  be  obtained  using 
only  the  coda  of  the  event,  individual  stations,  data  arriving 
through  individual  propagation  modes,  and  portions  of  the  data 
for  which  conventional  analysis  fails,  we  can  demonstrate  that 
this  new  procedure  extracts  and  properly  interprets  considerably 
more  depth  phase  information  then  previous  methods.  The  ability 
to  accomplish  this  will  enable  seismic  source  depth  determinations 
for  events  having  poorer  signal /noise  ratios  and/or  recorded  at 
fewer  stations  than  was  previously  possible. 

In  Figure  3.1  are  plotted  the  first  portion  of  the  seismograms 
recorded  from  the  Illinois  Event  of  11/9/68.  To  appreciate  what 
this  automated  analysis  procedure  must  accomplish  consider  Figure 
3.2.  Here  are  plotted  19  cepstrums  computed  from  consecutive  12.8 
second  data  samples  along  the  seismogram  recorded  at  WHZYK.  This 
shows  complex  and  changing  cepstrum  patterns,  most  of  which  contain 
depth  phase  information  which  must  be  constructively  utilized. 

This  is  primarily  accomplished  using  the  cepstrum  matched  filter 
technique  in  conjunction  with  differential  travel  time  information 
for  several  propagation  modes. 
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I The  Illinois  Event  was  first  analyzed  using  the  first  two 

minutes  of  data  recorded  at  all  six  stations.  The  output  of 
the  automated  analysis  is  plotted  in  Figure  3.3  as  cumulative 
CMF  output  versus  depth  in  km.  The  result  is  a very  clear 
dominant  peak  indicating  the  correct  source  depth  of  2G±2  km 
for  this  event.  However,  a more  impressive  result  is  seen  in 
Figure  3.4.  Here  the  same  analysis  was  performed  using  only  the 
second  minute  of  data  and  again  a clear  detection  of  the  event 
depth  is  obtained.  Through  the  omission  of  the  first  minute  of 
data  we  have  allowed  only  the  later  propagation  modes  (PP,  PPP, 
PcP)  to  contribute  to  this  result  and  have  clearly  demonstrated 
that  this  analysis  procedure  is  indeed  constructively  utilizing 
depth  phase  information  normally  ignored. 

Another  way  of  showing  that  depth  phase  information,  arriv- 
ing in  the  later  seismic  propagation  modes  is  contributing  to  the 
source  depth  estimate,  is  to  compute  the  CMF  output  versus  depth 
for  the  data  arriving  from  individual  propagation  modes.  In 
Figure  3.5a  through  3.5d  are  the  analysis  results  for  the  P,  PP , 
PPP  and  PcP  modes  individually  processed  using  two  minutes  of  data 
* from  each  of  six  stations.  In  each  of  these  plots  a clear  detec- 

tion of  the  correct  depth  is  obtained. 

The  ability  of  this  analysis  procedure  to  utilize  the 
differential  travel  times  information  of  later  propagation  modes 
has  another  important  advantage.  That  is,  a major  improvement  in 
the  ability  to  obtain  seismic  source  depths  from  data  recorded 
at  a single  station.  The  reason  is  that  normally  one  requires 
that  there  be  moveout  in  the  depth  phase  delay  times  over  a 
suite  of  stations  in  order  to  be  confident  that  the  depth  phase 
has  been  detected.  However,  in  this  new  analysis  procedure, 
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Cumulative  CMF  output  versus  depth  using  2 minutes  of  data 

FOR  ALL  SIX  STATIONS  FOR  ALL  MODES  COMBINED  (ILLINOIS  EVENT) 


DEPTH  (KM) 


tations  (Illinois  Even 
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VERSUS  DEPTH  USING  Z MINUTES  OF  D 
S FOR  PPP  MODE  ONLY  (ILLINOIS  Eve 


CNF  OUTPUT  VERSUS  DEPTH  USING  2 MINUTES  OF  DATA  FOR  ALL 
SIX  STATIONS  FOR  PcP  MODE  ONLY  (ILLINOIS  EVENT) 


i 


DEPTH  (KM) 


depth  phase  delay  time  moveout  occurs  along  the  coda,  as  different 
propagation  modes  contribute,  as  well  as  between  stations.  Thus, 
the  final  CMF  output  versus  depth  obtained  from  the  analysis  of 
a data  recorded  at  a given  station,  reflects  the  degree  to  which 
the  various  delay  times  are  in  agreement  with  a given  source 
depth . 

In  Figure  3.6a  through  3.6f  are  the  CMF  outputs,  using  one 
minute  of  data,  form  the  individual  analysis  of  data  recorded  at 
KN-UT,  MN-NV,  PCZBC,  WHZYK , NP-NT  and  FB-AK.  As  can  be  seen  from 
these  figures,  in  all  but  the  closest  station  KN-UT,  the  automated 
analysis  procedure  gives  very  good  indications  of  a '26  km  source 
depth  through  the  analysis  of  data  recorded  at  an  individual 
station.  Such  a capability  is  very  important  since  many  small 
magnitude  events  will  only  be  recorded  on  very  few  stations. 

From  the  analysis  of  the  Illinois  Event  using  this  new  auto- 
mated seismic  source  depth  determination  procedure,  all  indications 
are  that  the  concepts  and  techniques  used  in  the  analysis  arc 
valid  and  that  substantial  increases  in  the  percentage  of  deeper 
events  (>S  km)  for  which  source  depths  can  be  established  should 
be  realized. 

In  the  next  section  we  discuss  the  applicability  of  these 
techniques  to  shallow  events. 
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CMF  OUTPUT  VERSUS  DEPTH  USING  1 MINUTE  OF  DATA 
RECORDED  AT  STATION  KN-UT  (ILLINOIS  EVENT) 


t 


DEPTH  (KM) 


CMF  OUTPUT  VERSUS  DEPTH  USING  1 MINUTE  OF  DATA 
RECORDED  AT  STATION  fiN-NV  (ILLINOIS  EVENT) 
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DEPTH  (KM) 


CNF  OUTPUT  VERSUS  DEPTH  USING  1 MINUTE  OF  DATA 
RECORDED  AT  STATION  PGZBC  (ILLINOIS  EVENT) 


3-14 


f 


DEPTH  (KM) 


CKF  OUTPUT  VERSUS  DEPTH  USING  1 MINUTE  OF  DATA 
RECORDED  AT  STATION  WHZYK  (ILLINOIS  EVENT) 
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Figure  3.6d 


CNF  OUTPUT  VERSUS  DEPTH  USING  1 MINUTE  OF  DATA 
RECORDED  AT  STATION  NP-NT  (ILLINOIS  EVENT) 


t 


I 


I 
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f 


_i 


DEPTH  (KM) 


CHF  OUTPUT  VERSUS  DEPTH  USING  1 MINUTE  OF  DATA 
RECORDED  AT  STATION  FB-AK  (ILLINOIS  EVENT) 
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Figure  3.6f 


4.0  SHALLOW  SOURCF.  DEPTH  DETERMINATION  USING  CEPSTRUM  ANALYSIS 


In  order  to  apply  cepstrum  analysis  to  the  problem  of 
determining  delay  times  of  surface  reflected  seismic  energy  for 
shallow  source  depths  (<5  km),  one  faces  the  problem  of  short 
delay  times  ( <1  sec)  relative  to  the  available  signal  bandwidth 
(.5-3  Hz).  As  a consequence,  the  modulation  of  the  spectrum 
contains  less  than  a few  cycles  and  becomes  difficult  to  dis- 
tinquish  from  the  power  spectrum  of  seismic  arrivals  not  contain- 
ing surface  reflections. 

In  order  to  determine  the  effectiveness  of  the  automated 
analysis  and  to  shed  some  light  on  possible  improvements  and 
limits  of  cepstrum  analysis  for  these  shallow  depths,  we  analyzed 
both  real  and  synthesized  seismic  data.  This  included  the  analysis 
of  the  Boxcar  fvent,  known  to  be  at  a depth  of  1.12  km,  and  the 
Illinois  Event  having  a source  dentil  of  '26  km  for  the  nuroosc  of 
comparison.  This  analysis  was  done  using  the  automated  analysis 
procedure  modified  for  short  delay  times  by  the  removal  of  speefal 
tapering  and  the  cepstrum  matched  filter  (not  effective  in  its 
present  form  for  shallow  depths)  and  by  increasing  the  interpola- 
tion of  the  cepstrum. 


4 . 1 Analysis  of  the  Boxcar  Event 

Analysis  of  the  Boxcar  Event  was  done  using  data  recorded 
at  five  stations  having  epicenter  distances  of  15°,  17°,  21°, 
37°  and  38°.  Seven  12.8  second  consecutive  data  samples  from 
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each  station  recording  were  used  in  the  analysis.  The  first  50 
seconds  of  these  seismograms  are  plotted  in  Figure  4.1.1a. 

Figure  4.1.1b  is  a plot  of  the  averaged  amplitude 
spectrum  for  the  35  Boxcar  time  samples  over  the  frequency 
range  0-2.5  Hz. 

We  see  that  we  arc  dealing  with  a fairly  narrow  hand  signal 
which  in  effect  can  be  considered  to  have  a half  cycle  modulation 
in  the  amplitude  spectrum.  Figure  4.1.2  is  the  plot  of  the  cumu- 
lative CMF  output  versus  depth  in  km  resulting  from  the  modified 
automated  analysis  of  the  Boxcar  Event.  This  result  shows  little 
more  than  an  inflection  in  the  vicinity  of  the  known  depth  of 
1.12  km.  However,  if  one  computes  the  ceps t rum  using  Maximum 
Entropy  spectral  analysis  for  the  transform  of  the  amplitude 
spectrum,  the  results  of  Figure  4.1.3  show  a definite  peak  in 
the  1.1  km  range.  This  occurs  since  Maximum  Entrony 
spectral  analysis  is  known  to  be  more  effective  at  obtaining 
magnitudes  of  f-equency  components  when  less  than  a few  cycles 
arc  present.  One  variable  in  this  procedure  is  the  filter  length 
which  we  settled  on  125  out  of  a sample  containing  256  points. 

The  peak  at  '1.1  km  appears  to  be  a very  encouraging  result; 
however,  we  must  remember  that  the  shape  of  the  narrow  band  source 
spectrum  itself  can  introduce  structure  in  the  cepstrum  at  these 
short  delay  times.  We,  therefore,  performed  the  identical  analysis 
(using  Maximum  Entropy)  on  the  Illinois  Event  known  to  be  at  a 
greater  depth  and  having  a somewhat  similar  amplitude  spectrum 
as  seen  in  Figure  4.1.4.  The  result  of  this  analysis,  shown  in 
Figure  4.1.5,  shows  an  inflection  at  the  1.1  km  depth  range 
and  together  with  the  results  of  the  CMF  output  for  the  FF  mode 
seen  in  Figure  4.1.6,  one  sees  that  the  spectral  shape  can  indeed 
introduce  such  features  and  the  results  of  the  Boxcar  analysis 
must  be  interpreted  with  this  in  mind. 
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Figure  4.1.2 
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Figure  4.1.5 
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* Figures  4.1.7  and  1.1.8  show  the  results  of  the  analysis 

for  the  Boxcar  and  Illinois  F;vents  (not  using  Maximum  F.ntropy) 
with  the  depth  analysis  carried  out  to  60  km  and  the  mean  of  the 
amplitude  spectrum  removed.  One  notes  the  greater  structure  at 
depths  >10  km  from  the  Illinois  Event  compared  with  these  from 
the  Boxcar  Event.  One  sees  peaks  at  25  km  and  57  km  corresponding 
to  the  pP-P  and  sP-P  time  delays  for  the  Illinois  Event.  (The 
sP-P  delay  shows  up  since  the  CMF  was  not  used  in  the  run  and 
this  energy  was  not  folded  back  into  the  pP-P  peak.)  This  addi- 
tional structure  may  prove  to  be  a discriminate  for  natural  event  . 


4 . 2 Determination  of  the  Effectiveness  of  Ccpstrum  Analysis  for 

Shallow  Source  Depths  Using  Synthesized  Data 

In  order  to  determine  the  effectiveness  of  the  ccpstrum 
technique  in  determining  depths  for  shallow  events,  seismograms 
were  synthesized  to  simulate  events  having  surface  reflected 
delay  times  of  1.5,  .75,  .4  and  .25  seconds  corresponding  to 
depths  of  5,  2.5,  1.3  and  .83  km  for  epicenter  distances  of  '40°. 
This  was  accomplished  by  delaying  and  summing  Boxcar  seismograms 
to  themselves  (to  insure  realistic  signal  bandwidths  in  the  syn- 
thesized seismogram)  at  delays  corresponding  to  given  source 
depths . 

For  each  of  the  time  delays  used,  seismograms  were  generated 
for  five  stations  each  having  seven  coda  pieces  of  256  sample 
points.  The  automated  source  depth  analysis  procedure  was  used 
to  compute  an  average  weighted  ccpstrum  from  this  data.  A 
differential  ccpstrum  was  also  computed  by  scaling  and  subtract- 
ing from  these  ccpstrums  the  cepstrums  computed  from  the  seismo- 
grams not  containing  the  simulated  surface  reflection.  The 
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subtraction  technique  has  the  advantage  of  removing  those  ccpstrum 
features  which  are  generated  from  the  source  waveform  alone. 
Although  it  is  not  obvious  that  this  technique  can  be  used  with 
real  data  it  is  useful  in  determining  the  limitations  of  the 
cepstrum  technique  for  these  situations.  However,  with  additional 
information  such  as  surface  shots  and/or  use  of  the  theoretical 
relationships  between  waveforms  generated  at  different  depths 
and  yields,  such  a procedure  could  be  effective. 

Figures  4.2.1a  and  4.2.1b  show  tne  resulting  averaged  ccp- 
strum and  averaged  differential  cepstrum  for  a time  delay  of  1.5 
seconds.  The  peak  at  1.5  seconds  is  seen  in  both  cepstrums  but 
is  much  more  evident  in  the  differential  cepstrum  as  can  be 
expected  since  the  common  cepstrum  features  have  been  removed, 
leaving  only  the  feature  due  to  the  simulated  surface  reflection. 

, Figures  4.2.2a  and  4.2.2b  show  similar  plots  for  a time  delay  of 

.75  seconds.  Here  the  advantage  of  the  differential  cepstrum 
becomes  obvious  in  giving  a clear  peak  at  .75  seconds  versus  only 
an  inflection  i:i  the  cepstrum  computation.  Figures  4.2.3a  and 
( 4.2.3b  show  a case  in  which  the  time  lag  has  been  reduced  to  .4 

seconds,  where  :here  is  no  indication  of  the  time  lag  in  the 
averaged  cepstrum  but  is  again  very  clear  in  the  differential 
cepstrum.  Finally,  we  used  a delay  of  .25  seconds  and  plots  of 
i Figures  4.2.4a  and  4.2.4b  indicate  that  the  differential  cepstrum 

is  still  capable  of  detecting  this  time  lag  but  the  signal/noise 
ratio  has  greatly  deteriorated;  primarily  a consequence  of  the 
short  delay  and  bandlimited  signal.  The  analyses  of  these  syn- 
i thesized  seismograms  were  performed  to  indicate  what  one  can 

expect  from  cepstral  analysis  for  shallow  depth  phase  detection 
for  realistic  signal  bandwidths.  A result  of  this  research  is 
that  one  needs  to  develop  techniques  which  allow  the  estimation 
» of  cepstrum  features  originating  from  the  source  waveform  not 

containing  the  surface  reflection  in  order  to  detect  delay  times 
of  <1  second. 
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Figure  4.2.1a 
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Figure  4.2.2a 
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Figure  4.2.2b 
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Figure  4.2.3a 
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Figure  4.2.4a 


4-20 


5.0  CONCLUSION 


During  this  contract,  an  automated  tcleseismic  source  depth 
determination  analysis  has  been  developed  for  events  having  source 
depths  in  the  50-200  km  range,  and  was  applied  to  seismic  event 
data.  It  was  demonstrated  that  through  this  analysis,  source 
depth  information  contained  throughout  the  seismogram  is  being 
constructively  utilized  to  obtain  source  depth  estimates. 

It  was  also  demonstrated  that  this  analysis  gives  correct  source 
depths  from  data  recorded  at  individual  stations  for  the 
Illinois  Event.  These  capabilities  should  allow  one  to  obtain 
source  depths  for  a considerably  higher  percentage  of  seismic 
events . 

Research  was  conducted  to  establish  possible  modifications 
to  the  analysis  which  would  enhance  its  effectiveness  for  events 
having  source  depths  less  than  5 km.  A result  of  this  research 
was  that  in  order  to  use  the  cepstrum  analysis  approach  for 
shallow  events,  techniques  to  determine  cepstrum  features  result 
ing  from  the  direct  seismic  arrival  (those  not  containing  the 
surface  reflected  arrival)  need  to  be  developed. 

To  establish  the  general  applicability  and  conditional 
limitations  of  the  source  depth  analysis,  a large  set  of  events 
must  be  analyzed.  In  addition,  techniques  to  establish  the 
significance  and  estimated  error  in  these  depth  determinations 
should  be  included  in  the  analysis  procedure. 
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APPENDIX  A 

COMPUTATION  OF  CEPSTRUM 


» 


The  following  steps  are  used  to  calculate  the  Stochastic 
Cepstrum  for  a seismogram  consisting  of  NS  consecutive  portions 
each  having  N amplitude  values: 

• Select  N sampled  (sample  rate  = 20  pts/sec) 

X.  h 

amplitude  values  from  the  portion  of  a 
seismegram 

XM(n) , n=0 , 1 , . . . , N-l 
and  add  N zeroes  to  interpolate  the  spectrum 

• Calculate  the  amplitude  spectrum  for  positive 
frequencies  using 


2N-1 


A„0)  - 


_i_ 

2N 


^ XM(n)e-2*ini/2N 


; j-0,1,...,  N-l 


n=0 


1/2 

with  the  Nyquist  frequency  = 10  Hz  and  i = (-1)  ' . 

• Retain  only  the  lower  quarter  frequencies  of 
the  amplitude  spectrum  since  there  is  very 
little  energy  at  frequencies  above  2.5  llz 
for  natural  events.  This  leaves  the  array 


A- 1 


(AM(j),  j = 0,  1,  ....  N/4-1) 


• Remove  the  mean  and  apply  a cosinusoidal  taper 
to  the  first  10%  and  last  20 % of  the  array 
giving  the  modified  array 

(Am  (j)  , j - 0,  1,  N/4-1) 

The  201  taper  on  the  higher  frequencies  was  used 
to  de-emphasize  the  higher  frequencies 
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• Add  N/4  zeroes  to  interpolate  the  cepstrum  giving 
the  array 

(Am  (j)  , j = 0,  1 N/2-1) 

At  this  stage  one  would  take  the  log  of  this 
A'  array  to  obtain  the  log  cepstrum;  but  by  not 
using  the  log  better  results  were  obtained. 


• Calculate  the  Fourier  Transform  of  the  AM  array 
using 


N/2-1 

FM<k)  ■ m Z am  (i)e 

j-0 


..  -2irijk/(N/2) 


A-  2 


I 


I 


where  (F^(k),  k = 0,  1 N/4-1)  are 

complex  numbers  representing  the  positive  fre- 
quency spectrum  of  A^.  One  can  now  obtain  the 
amplitude  and  phase  of  each  copstrum  point  k. 

The  array 

( | F^j  ( k ) | , k = 0,  1,  N/4-1)  is  what  we 

is  referred  to  as  the  cepstrum  amplitude  and  the 
complex  numbers  F^(k)  are  referred  to  as  cepstrum 
phasors  . 

• To  calculate  the  stochastic  cepstrum,  one 
then  calculates 

|FM(k)|  = Max  ( |PM(j)|,  j = k-A/2,  k+A/2) 

for  each  section  M and  sums  these  over  the  number 
of  sections  (NS)  used  from  a given  seismogram 
(A  is  the  stochastic  window  width).  This  gives 


NS 

CCk)  = ^ . 

M=1 

where  C(k)  is  the  stochastic  cepstrum  and 
WM  is  a weighting  factor,  chosen  such  that  the 
mean  amplitude  of  each  |F^(k)|  array  is  equal. 
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APPENDIX  B 

CEPSTRUM  MATCHED  FILTER  (CMP) 


If  one  assumes  that  the  cepstral  peaks  appearing  at  the 
pP-P,  sP-P  and  sP-pP  time  delays  dominate  the  expected  cepstrum 
pattern  for  an  event,  then  one  can  formulate  the  CMP  algorithm 
in  the  following  way.  Consider  the  cepstrum  to  consist  o<"  N 
amplitude  values  CP(n)  for  the  time  delays  t =(n-l)*At,  where 
n = 1,2,3...N.  Then  the  synthetic  cepstrum  (CS)  expected  for  an 
event  at  a given  depth  and  distance  corresponding  to  a pP-P  delay 
1 time  of  t,  con  he  represented  by: 

CS (n , t , a)  = Q(n,(a-l)x)  + Q (n , t ) + Q(n,ax) 

f 

witha=  Tsp.p/Tpp_p,  and  is  the  ratio  of  delay  times  for  sP-P 
and  pP-P.  Also  C is  defined  by 

Q(n,x)  = 1 ; if  t - At<tn<x+At  , 

= 0 ; otherwise  . 

Thus,  CS  is  represented  by  three  unit  amplitude  peaks  located  at 
the  sP-pP,  pP-P,  and  sP-P  delay  times,  each  peak  consisting  of 
three  adjacent  delay  time  points. 

The  CMP  output  at  time  delay  t is  defined  to  be  the  maximum 
zero  lag  correlation  of  the  computed  cepstrum  (CP)  with  the 
synthesized  cepstrum  (CS)  for  a source  depth  having  a pP-P  delay 
time  x computed  over  a range  of  ratios  a.  This  can  be  written 

as 
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CMF(t) 


i f CP (t ) > . 7 • CP  (a t ) 
and  CP(t)>.7‘CP((o-1Jt) 


Max 

a j <_  a <_  c<2 


)'  C P (n)  • CS  (n  , t , a ) ' 
n=  1 


= CP(t)  ; otherwise 


Here,  a.  and  <*7  are  set  from  the  expected  range  of  values 
for  a for  a given  range  of  possible  depths.  Values  of  = 1.25 
and  « 2 = 1.55  were  used  for  the  event  distances  and  depths  encoun- 
tered in  this  research.  The  constraints  imposed  require  that  the 
amplitude  of  the  cepstrum  peak  at  t must  he  at  least  .7  times 
the  amplitude  of  the  cepstrum  peaks  at  both  ai  and  (cx-1)t  • 

This  eliminates  detection  of  cepstral  patterns  for  cases  in 
which  the  strength  of  the  pP  arrival  is  considerably  less  than 
that  of  the  sP  arrival.  Thus,  a cepstrum  peak  at  t resulting 
from  a pP  arrival  for  an  event  will  not  contribute  to  the  CMF 
. output  for  an  apparent  event  having  an  expected  sP-P  delay  time 

of  t,  without  the  significant  presence  of  a pP  arrival  for  this 
apparent  event. 

As  an  example  of  the  application  of  this  technique,  consider 
Figure  1.  At  the  top  of  the  figure  is  the  cepstrum  calculated 
from  an  event  of  known  depth  having  both  the  pP  and  sP  depth 
phases  clearly  identifiable.  This  cepstrum  pattern  is  dominated 
by  three  peaks  corresponding  to  the  sP-pP,  pP-P,  and  sP-P  delay 
times.  Upon  interpreting  this  cepstrum  for  an  unknown  event  depth, 
one  secs  that  each  of  the  three  peaks  has  the  possibility  of 
corresponding  to  the  pP  phase.  The  CMF  output  is  plotted  at  the 
bottom  of  Figure  1 and  indicates  a much  stronger  emphasis  on  the 
peak  corresponding  to  the  correct  pP-P  delay  time  for  this  event. 
This  result  reflects  the  fact  that  this  cepstrum  pattern,  primarily 
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consisting  of  three  peaks,  is  in  strong  agreement  with  the  pattern 
expected  for  a single  event  having  this  pP-P  delay  time.  The 
lesser  probability  that  the  relative  location  of  these  three  peaks 
was  coincidental,  and  that  one  of  the  other  two  peaks  actually 
corresponded  to  the  correct  pP-P  delay  time,  is  indicated  by  the 
presence  of  CMF  output  peaks  at  those  delay  times  having  reduced 
ampl i tudes . 


» 


» 
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Figure  1.  Ceps t rum  ami  Ccpst rum  Matched  Filter  output 

computed  from  data  recorded  at  several  stations 
for  an  event  with  well  defined  depth  phases. 
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